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Abstract 

We report polarimetry results of a merging cluster of galaxies Abell 2256 with Karl G. Jansky 
Very Large Array (JVLA). We performed new observations with JVLA at S-band (2051- 
3947 MHz) and X-band (8051-9947 MHz) in the C array configuration, and detected significant 
polarized emissions from the radio relic, Source A, and Source B in this cluster. We calculated 
the total magnetic field strengths toward the radio relic using revised equipartition formula, 
which is 1.8-5.0 pG. With dispersions of Faraday rotation measure, magnetic-field strengths 
toward Sources A and B are estimated to be 0.63-1.26 /jG and 0.11-0.21 pG, respectively. An 
extremely high degree of linear polarization, as high as ~ 35 %, about a half of the maximum 
polarization, was detected toward the radio relic, which indicates highly ordered magnetic lines 
of force over the beam sizes (~ 52 kpc).The fractional polarization of the radio relic decreases 
from ~ 35 % to ~ 20 % around 3 GHz as the frequency decreases and is nearly constant 
between 1.37 and 3 GHz. Both analyses with depolarization models and Faraday tomography 
suggest multiple depolarization components toward the radio relic and imply the existence of 
turbulent magnetic fields. 

Key words: galaxies: clusters: individual: (Abell 2256) — galaxies: clusters: intracluster medium — 
magnetic fields — polarization 


1 Introduction 

Collision of galaxy clusters is one of the most energetic events 
with kinetic energy on the order of ~ 10 64 ergs in the 

Universe (Buote 2001). Shock waves and turbulence induced 

* Mail to: k5148776@kadai.jp 


by the collision can convert a huge kinetic energy of clusters 
into thermal/non-thermal energies of the intracluster medium 
(ICM). Cosmic-rays injected into the ICM by AGN activi¬ 
ties, star formations of galaxies, and structure formation shocks 
(Brunetti et al. 2001) can be re-accelerated by the shocks 
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(Takizawa & Naito 2000; Ryu et al. 2008; Vazza et al. 2009) 
and turbulence (Brunetti et al. 2001; Petrosian 2001; Ohno et 
al. 2002; Fujita et al. 2003; Brunetti et al. 2004; Cassano & 
Brunetti 2005; Xu et al. 2009; Xu et al. 2010; Feretti et al. 
2012; Donnert et al. 2013). The correlation between X-ray 
luminosity of the ICM and the power of diffuse radio emis¬ 
sion from cosmic-rays is known for radio halos and relics (e.g., 
Feretti et al. 2012). It suggests the relationship between the 
cluster size and the magnitude of particle acceleration in sense 
that larger clusters can produce more powerful radio emission. 
However, the nature and evolution of the ICM and intergalactic 
magnetic field (IGMF), which determine the efficiency and the 
radio emission mechanisms, are poorly understood. 

Turbulence is thought to play an important role in the evo¬ 
lution of the IGMF. It has been suggested that turbulence dy¬ 
namo can amplify the IGMF in a cosmological time (Ryu et al. 
2008; Cho & Ryu 2009). Actually, the Kolmogorov index in 
the power spectrum of magnetic fields has been reported (e.g. 
Abell 2362, Guidetti et al. 2008), indicating the existence of tur¬ 
bulence and amplification of the IGMF by turbulence in galaxy 
clusters. 

One of the useful techniques to investigate turbulent mag¬ 
netic fields is the depolarization, in which observed polarized 
intensity gets weaker than that at the origin arising from sev¬ 
eral mechanisms. Particularly, beam depolarizations of inter¬ 
nal Faraday dispersion (IFD) and external Faraday dispersion 
(EFD) becomes significant effect, if structures of Faraday rota¬ 
tion measure (RM) inside and outside a polarized radio emis¬ 
sion source, respectively, are not uniform within an observing 
beam, e.g. if small eddy-sized turbulent magnetic fields exist. 
Both IFD and EFD depend on the dispersion of RM within the 
beam. Burn (1966) analytically investigated the dependencies 
called Burn’s law, and Arshakian & Beck (2011) investigated 
the optimum frequency range for this technique. 

Abell 2256 is known as a merging cluster of galaxies in 
which we expect turbulence of the ICM and turbulent magnetic 
fields in the cluster. In this paper, we report results of linear 
polarimetry of the central part of Abell 2256 with the Karl G. 
Jansky Very Large Array (JVLA) at S-band (2051-3947 MHz) 
and X-band (8051-9947 GHz) in the C array configuration. We 
obtained Stokes I, Q, and U images in order to measure the 
total intensity, fractional polarization, and RM map. Using the 
obtained maps, we investigated IGMF structures in Abell 2256 
by means of depolarization. The layout of this paper is as fol¬ 
lows. In Section 2, we introduce Abell 2256. In Section 3, we 
describe the observations and data reductions in the JVLA. In 
Section 4, we present the results, which include the total inten¬ 
sity, fractional polarization and RM. In Section 5, we discuss 
the magnetic field strengths toward the radio relic, Source A, 
and Source B, and discuss the fractional polarization of the ra¬ 
dio relic. In Section 6, we summarize our conclusions. 


Throughout this paper, we assume the following cosmolog¬ 
ical parameters: Ho = 70.5 km s -1 Mpc -1 , D m o = 0.27, and 
Dao = 0.73. The angular size of 1' corresponds to ~ 67 kpc 
at the redshift of Abell 2256, 2 = 0.0581, corresponding to a 
distance of D = 247 Mpc. 

2 Cluster of galaxies Abell 2256 

Abell 2256 is a nearby (redshift 2 = 0.0581) cluster of galax¬ 
ies whose X-ray center is located at (RA, Dec) = (17 h 04 m 2?3, 
+78°37'55."2) in the J2000 epoch (Ebeling et al. 1998). 
Berrington et al. (2002) investigated Abell 2256 with opti¬ 
cal observations and found substructures of member galaxies 
with the peak radial velocity difference of ~ 2000 kmT 1 . 
Substructures of the ICM are also known in X-ray observations 
(Briel et al. 1991; Briel & Henry 1994; Sun et al. 2002). Using 
the X-ray satellite Suzaku, Tamura et al. (2011) estimated radial 
velocity difference of ~ 1500 km s' 1 in gas bulk motions of 
the substructures. There are two distinct ICM components with 
temperatures ~ 7 keV and ~ 4.5 keV (Sun et al. 2002). These 
results suggest that Abell 2256 is a merging galaxy cluster. 

Radio observations have discovered a radio relic and a 
halo in the central part of Abell 2256 (Bridle & Fomalont 
1976; Bridle et al. 1979; Rottgering et al. 1994; Miller et al. 
2003; Clarke & Ensslin 2006; Brentjens 2008; van Weeren 
et al. 2009; Kale & Dwarakanath 2010; van Weeren et al. 
2012; Owen et al. 2014; Trasatti et al. 2015). The radio relic 
located in the north-west of the cluster is ~ 440 kpc away 
from the X-ray center. The radio relic covers an area of 16.'9 x 
7'.8 (1125 kpc x 520 kpc) (Clarke & Ensslin 2006). Previous 
observations revealed that the radio relic includes filamentary 
structures (Clarke & Ensslin 2006; Brentjens 2008; Owen et al. 
2014). The radio halo is located in the central region of the clus¬ 
ter. Clarke & Ensslin (2006) measured that the total flux of the 
radio halo is approximately 103 mjy at 1369 MHz. 

There are also several radio sources in Abell 2256. Miller et 
al. (2003) identified the radio sources associated with member 
galaxies of Abell 2256. Each radio source is labeled with alpha¬ 
bets (Bridle & Fomalont 1976; Bridle et al. 1979; Rottgering et 
al. 1994), and Sources A, B, and C are the remarkable bright 
sources. Sources A and B are linearly polarized sources which 
are suitable for our measuring RMs. Source C is known as a 
head-tail galaxy which has a narrow straight tail extending for 
at least 480 kpc at 1.4 GHz across the radio relic (Rottgering et 
al. 1994). 

3 Observations and Data reductions 

3.1 Radio observations 

We carried out new observations of Abell 2256 using JVLA 
at S-band (2051-3947 MHz) and X-band (8051-9947 MHz) in 



Publications of the Astronomical Society of Japan, (2014), Vol. 00, No. 0 


3 


Table 1. Details of the VLA & JVLA observations of Abell 2256. 



Frequency* 

(MHz) 

Bandwidth* 

(MHz) 

Configuration * 

Date 

Time* 

(h) 

Project* 

VLA 

1369/1417 

25/25 

D 

1999-Apr-2 8 

5.9, 5.9 

AC0522 


1513/1703 

12.5/25 

D 

1999-Apr-29 

3.5, 5.5 


VLA 

1369/1417 

25/25 

C 

2000-May-29 

2.5, 2.5 

AC0545 


1513/1703 

12.5/12.5 

C 

2000-May-29 

3.6, 3.6 



1369/1417 

25/25 

C 

2000-Jun-18 

2.5, 2.5 



1513/1703 

12.5/25 

C 

2000-Jun-18 

4.1, 3.5 


JVLA 

S-band 16 windows^ 

128 

C 

2013-Aug-25 

1.2 

13A-131 (this work) 





2013-Aug-26 

1.2 






2013-Aug-29 

1.2 


JVLA 

X-band 16 windows 1 *- 

128 

C 

2013-Aug-18 

1.3 

13A-131 (this work) 





2013-Aug-19 

1.3 



* Column 2: observing frequency; Column 3: observing bandwidth; Column 4: array configuration; Column 5: dates of observation; Column 
6: time on source; Column 7: NRAO project code. 

1 2051/2179/2307/2435/2563/2691/2819/2947/3051/3179/3307/3435/3563/3691/3819/3947. 

* 8051/8179/8307/8435/8563/8691/8819/8947/9051/9179/9307/9435/9563/9691/9819/9947. 


the C array configuration on 18-30 August 2013. Each band 
was separated into 16 spectral windows and each window had 
a bandwidth of 128 MHz (Table 1). The pointing center was 
Source A (17 h 03 m 3L 5 9,+78 o 37 , 44"4), since our observations 
primarily aimed at measuring RMs toward Sources A and B 
in the central part of Abell 2256. We observed 3C286 and 
1803+784 as a flux and polarization calibrator, and a gain and 
phase calibrator, respectively. 

3.2 Data reductions 

Data were reduced by the following procedures. Using National 
Radio Astronomy Observatory (NRAO) Common Astronomy 
Software Applications (CASA), we executed VLA calibra¬ 
tion pipeline and task EXPORTUVFITS to convert JVLA’s 
measurement sets into FITS so as to allow us to calibrate 
the data with NRAO Astronomical Image Processing System 
(AIPS). In AIPS, the data in each spectral window was aver¬ 
aged in frequency domain using task AVSPC. Radio Frequency 
Interferences (RFls) and spurious signals were flagged using 
task CLIP and TVFLAG. We calibrated polarization using task 
PCAL and RLDIF after we created calibration tables. The 
data separated into observed date are concatenated using task 
DBCON. 

The data at 2179, 2307, 3691, 3819, and 3947 MHz in S- 
band are affected by RFIs. There were satellite downlink and 
digital audio radio service in 2180-2290 MHz and 3700M-200 
MHz, respectively 1 . Therefore, we removed the data in these 
spectral windows from our analysis. 

In addition to the observed data, we also utilized the archival 
data which were observed with the VLA in the C and D array 
configurations at L-band (1369, 1417, 1512, 1703 MHz). We 
calibrated these data by the same procedures described above. 

1 https://science.nrao.edu/facilities/vla/docs/manuals/oss/performance/rfi 


We created Stoke I, Q and U images using task IMAGR 
with suitable tapers. In order to detect the radio relic and re¬ 
solve the individual polarized sources in the cluster, we made 
images of 47" resolution for L and S bands, and 15" 1 resolution 
for S and X bands (Table 2). Note that in the following analy¬ 
ses except Section 4.1, we used the images of 47" and 15."1 
resolution for analyzing the radio relic and the individual radio 
sources, respectively. We excluded the data of X-band from the 
images of 47" resolution, since the radio relic is outside the field 
of view. The images was convolved with a Gaussian beam using 
task CONVL. The number of pixels was 87 x 87 for images of 
47" resolution, and 136 x 136 for images of 15" 1. Each pixel 
size corresponds to a half size of each beam size. 


4 Results 

4.1 Radio images 

Figure 1 shows the total intensity map of Abell 2256 at 2051 
MHz. We detected several radio sources and the radio relic 
above 3ai significance. All of them are known sources in previ¬ 
ous works with different array configurations and/or frequency 
bands (e.g. Miller et al. 2003; Clarke & Ensslin 2006; Owen 
et al. 2014). In all of created images in S-band except the data 
affected by RFIs, we detected the radio relic, where an area of 
the relic above 3cti significance at 3563 MHz was smaller than 
that at 2051 MHz by a factor of ~ 0.37. 

Figure 2 shows the total intensity map at 8051 MHz. The 
radio relic is outside the field of view. In all of created images in 
X-band, we detected Sources A, B, and C. These radio sources 
are identified as radio galaxies (Miller et al. 2003). 
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Table 2. Image qualities of total intensity and polarization at L, S, and X bands. 


Frequency* 

(MHz) 

Beam* 

("x") 

OI* 

(mjy beam -1 ) 

^Q* 

(mjy beam -1 ) 

(mjy beam -1 ) 

VLA 1369 

47x47 

0.163 

0.028 

0.022 

1417 

47x47 

0.152 

0.031 

0.020 

1513 

47x47 

0.183 

0.041 

0.028 

1703 

47x47 

0.259 

0.046 

0.095 

JVLA S-band 11 windows^ 

47x47 

0.159 

0.029 

0.029 

S-band 11 windows^ 

15.1x15.1 

0.151 

0.013 

0.014 

X-band 16 windows $ 

15.1x15.1 

0.053 

0.027 

0.028 


* Column 2: observing frequency; Column 3: beam size; Columns 4,5,6: RMS noise of the Stokes I, Q, and U. We show the 
averaged RMS noise in JVLA S-band 11 windows and X-band 16 windows. 

+ 2051/2435/2563/2691/2819/2947/3051/3179/3307/3435/3563. 

1 8051/8179/8307/8435/8563/8691/8819/8947/9051/9179/9307/9435/9563/9691/9819/9947. 



Fig. 1. Total intensity map of Abell 2256 at 2051 MHz in the JVLA C ar¬ 
ray configuration. Contours are drawn at (-3,3,6,12,24,48,96,192) x 
140 n Jy beam -1 (9.702 X 10 -23 W m -2 Hz -1 sr -1 ). The beam size 
of the image is shown at bottom-left and is 25."8x21."0. Radio sources 
are labeled following the references (Bridle & Fomalont 1976; Bridle et al. 
1979; Rottgering et al. 1994). The fractional polarization of the radio relic 
was measured in the dashed frame region (Section 4.3). 


4.2 Total intensity 


Fig. 2. Total intensity map of Sources A, B, and C at 8051 
MHz in the JVLA C array configuration. Contour levels are 
drawn at (-3, 3, 6, 12, 24, 48, 96, 192) X 41.2 /xjybeam -1 
(9.424 x lo -23 W m -2 Hz -1 sr -1 ) . The beam size of the image 
is shown at bottom-left and is 14."4x11."4. 

~ 13 %. Note that the total 2051 MHz flux of 286 mjy is not 
subtracted the flux from the tail of Source C, so that the total 
flux of the relic is smaller than 286 mjy. 


The total 2051 MHz flux of the radio relic is estimated to be 
286 mjy based on the image. For calculating the total 2051 
MHz flux, we integrated the pixels where the flux density is 
above 3 <j/ significance on the radio relic in the image of 47" 
resolution. If we consider the total 1369 MHz flux from the 
radio relic of 462 mjy which is subtracted the flux of the point 
sources (Clarke & Ensslin 2006), and adopt the spectral index 
(S oc v a ) of a = —0.81 (van Weeren et al. 2012), the total 2051 
MHz flux should be 330 mjy. The measured total 2051 MHz 
flux of 286 mjy is thus smaller than the expected total flux by 


We also estimated the upper limit of the total 2051 MHz 
flux of the radio halo, which is ~ 49 mjy, assuming that 
the radius of the radio halo emission is ~ 6(1 (Clarke & 
Ensslin 2006) and the upper limit of the flux density is 
606 /rJy beam -1 (10.296 x 10 - 23 W m -2 Hz -1 sr -1 ) which 
corresponds to 3or significance at 2051 MHz. If we adopt the 
total 1369 MHz flux of the radio halo of 103 mjy (Clarke & 
Ensslin 2006) and the spectral index of a = —1.1 (van Weeren 
et al. 2012), the total 2051 MHz flux should be ~ 66 mjy. Thus, 
the estimated upper limit of the total flux of 49 mjy is smaller 
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Fig. 3. An example of the spectral energy distribution (SED). Open squares 
show the SED for a pixel (l7 h 02 m 51?9, +78°42'26."7) in the radio relic. 
The solid line represents an extrapolation of the spectrum from the surface 
brightness between 1369 MHz to 2051 MHz assuming the spectral index 
a = -0.81. The dashed line represents a measured spectral index of a = 
-1.98 from the surface brightness between 1369 MHz to 3563 MHz. 


than the expected total flux by ~ 26 % 

Figure 3 shows an example of the spectral energy distribu¬ 
tion (SED) at a point (17 h 02 m 51?9,+78°42'26"7) in the radio 
relic, where the SED is made from the images of 47" resolu¬ 
tion. We obtained a spectral index of a ~ —1.98 from the sur¬ 
face brightness between 1369 MHz to 3563 MHz (the dashed 
line in Figure 3). We found that the observed flux density of 
the radio relic at the frequency above 2 GHz is smaller than the 
extrapolation of the flux density from the results in 1369-2051 
MHz with the spectral index a = —0.81 (the solid line in Figure 

3). 


Fig. 4. The fractional polarization spectra of the radio relic (open squares). 
The filled square at the bottom-left is a observational result for the radio relic 
with the WSRT at 350 MHz (Brentjens 2008). The lines show the depolariza¬ 
tion models defined in equation (9) with p 0 =0.73, cr RM =80 rad m _2 (solid 
line) and p 0 — 0.2, o- RM = 6 rad m _2 (dashed line). 


0 . 
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Frequency [GHz] 


Fig. 5. The fractional polarization spectra of the Source A (open circles), 
Source B (open triangles), and Source C (filled inverted triangles). 


A possible cause of such a decline is that we are missing 
the flux. In interferometry, we miss the flux from the diffuse 
emission which has a scale larger than the largest angular scale 
(LAS) of the interferometer. Actually, the scale of the ma¬ 
jor axis of the radio relic and halo is 1014" and 732" on the 
sky (Clarke & Ensslin 2006), while the LASs at 1.5 GHz and 
3.0 GHz in the JVLA C array configuration were 970" and 
490", respectively. This possibility is also supported by single¬ 
dish radio observations using the Green Bank (Owen 1975) and 
Effelsberg (Haslam et al. 1978) telescopes, which yielded the 
total flux of 570 and 666 mjy in the entire area of Abell 2256 
at 2695 MHz, respectively, in agreement with the spectral index 
of a = —0.72 (Brentjens 2008). 

Another possible cause of the decline would be a cutoff of 
cosmic-ray electrons at high energies. But since the effect of 
missing flux could be significant, we could not argue the possi¬ 
bility of the energy cutoff. 


4.3 Fractional polarization 


We detected significant polarized emission from the radio relic, 
Source A, and Source B at S-band. At X-band, we detected 
the significant polarized emission only from Source A. With the 
Stokes I, Q, and U, the fractional polarization p is given by 


y /Q 2 + U 2 


(1) 


The fractional polarization are created from the images of 47" 
and 15!' 1 resolution to analyze the radio relic and the individ¬ 
ual polarized sources. We calculated the fractional polarization 
only in the pixels where the flux densities of Stokes I, Q, and 
U are all above 3 o significance. 

Figure 4 shows the fractional polarization spectra of the ra¬ 
dio relic (open squares), and Figure 5 shows the fractional polar¬ 
ization spectra of Source A (open circles), Source B (open trian¬ 
gles), and Source C (filled inverted triangles). Each data point 
represents a spatial average for the pixels within each emitting 
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region. For the radio relic, we choose the region where the po¬ 
larized emission is detected at 3563 MFlz (the dashed frame re¬ 
gion in Figure 1). We only plotted the data points which satisfy 
that the fractional polarization was obtained with at least 3 pix¬ 
els within each emitting region in each frequency. The error bar 
indicates the standard deviation of the fractional polarization for 
the pixels. We found that the fractional polarization of the ra¬ 
dio relic decreases from ~ 35 % to ~ 20 % as the frequency 
decreases from ~ 3.5 GHz to ~ 3 GHz. The fractional polariza¬ 
tion is then nearly constant between 1.3-3 GHz. We plotted the 
fractional polarization of the brightest part of the radio relic ob¬ 
served at 350 MHz with Westerbork Synthesis Radio Telescope 
(WSRT) (Brentjens 2008). It indicates that the fractional polar¬ 
ization is less than 1 %. Therefore, the fractional polarization 
in the field toward the relic varies twice, at ~ 3 GHz and around 
~ 0.4-1.3 GHz, and has step-like variations. 

We also see a change of the fractional polarization be¬ 
tween 8051 MHz to 9947 MHz for Sources B and C. This is 
due to an artifact because sensitivity is not enough at X-band. 
For instance, the polarized intensity (0.152 mjy beam -1 ) is 
smaller than 3cti significance (0.159 mjy beam -1 ) for a pixel 
(17 h 03 m 07 s 8, +78°36'17."8) in the Source B, and the polarized 
intensity (0.126 mjy beam -1 ) is smaller than 3or significance 
(0.145 mjy beam -1 ) for a pixel (17 h 03 m 28 s l, +78°39'56"2) 
in the Source C at 8051MHz. We consider that weak polarized 
emission exist at X-band but we did not detect the accurate po¬ 
larized intensity due to low sensitivity. 

4.4 Faraday Rotation Measure 

Faraday rotation measure (RM) is defined from the Faraday ro¬ 
tation of the linear polarization as 

<j> = ^> o + RMA 2 , (2) 

where tj> is the observed polarization angle in radians, and </>o is 
the intrinsic polarization angle in radians. RM is given by 

RM « 812 J n e B n dl, (3) 

in rad m -2 , where n e is the thermal electron density in cm -3 , 
B is the magnetic fields parallel to the line of sight in pG, and l 
is the path length in kpc. 

The RM map was created according to the linear relation be¬ 
tween cj> and A 2 in equation (2). We used the polarization angle 
images of 47" and 15" 1 resolution to analyze the radio relic and 
the individual polarized sources. We calculated RM only in the 
pixels which satisfy the following conditions: the flux density 
of the Stokes I, Q, and U are all above 3<r significances, and the 
pixels satisfying the first condition are available from at least 4 
frequencies. The RM map of Abell 2256 is shown in Figure 6 
and 7. 

To make sure that our RMs based on a linear-fit between (f> 



Fig. 6. The RM distribution of the Abell 2256 with the intrinsic B-vector. 
Also overlaid as contours is the total intensity map at 2051 MHz in the 
JVLA C array configuration drawn at (-3, 3, 6, 12, 24, 48, 96, 192) x 
201.9 fijy beam -1 (3.432 X 10 -23 W m -2 Hz -1 sr -1 ), convolved with 
a circular Gaussian beam with a FWHM of 47". The numbers represent the 
positions chosen in Figure 8. 



Fig. 7. The RM distribution of the Sources A, B, and C. Also overlaid as 
contours is the total intensity map at 8051 MHz in the JVLA C array con¬ 
figuration drawn at (—3,3,6,12,24,48,96,192) X 48.4 /rjy beam -1 
(7.976 x 10 -23 W m -2 Hz -1 sr -1 ), convolved with a circular Gaussian 
beam with a FWHM of 15."1. 
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Fig. 8. Sample plots of the polarization angle <j> against A 2 for different posi¬ 
tions in Abell 2256. Each position is shown in Figure 6 and 7. 


Table 3. The average and standard deviation of RM. 


Target 

(RM)* 
rad m —2 

CTRM* 
rad m —2 

reference 

Relic 

-44 

7 

Clarke & Ensslin (2006) 

Relic 

-34.5 

6.2 

this work 

Source A 

-24.9 

65.5 

this work 

Source B 

-34.1 

10.5 

this work 


* (RM) and ctrm are the average and standard deviation of RM. respectively. 


and A 2 are reasonable, we examined the <p-\ 2 relations toward 
the radio relic, Source A, and Source B. Figure 8 shows some 
examples for the positions inside them (see Figure 6 and 7). We 
confirmed that the linear relation is roughly satisfied for the ra¬ 
dio relic, Source A, and Source B. We also detected the RMs 
from the Source C but we did not use these RMs, since the po¬ 
larized emission from Source C is not significant at X-band due 
to low sensitivity (see Section 4.3). 

We calculated the average, (RM), and the standard devia¬ 
tion, ctrm, of RM for the radio relic. Source A, and Source B. 
The results are listed in Table 3. We also show the results for the 
radio relic reported by Clarke & Ensslin (2006). Our results for 
the radio relic and Source B are broadly consistent with the pre¬ 
vious estimations for the radio relic. On the other hand. Source 
A has substantially smaller (RM) and much larger ctrm com¬ 
pared to the other positions. Murgia et al. (2004) reported that 
the simulated |(RM)|/ctrm ratio depends only on the magnetic 
field power spectrum slope, and it has a considerable scatter. We 
consider that the smaller (RM) and larger ctrm could be related 
to the magnetic field fluctuations in the cluster. 


5 Discussion 

5.1 Magnetic Field Strengths in the Radio Relic 

We calculated the magnetic field strengths of the radio relic in 
Abell 2256 using the revised equipartition formula from Beck & 
Krause (2005). The total equipartition magnetic field strengths 


B t is given by 

4tt( 1 - 2a)(_K'o + l)BEl~ 2a {v /2ci)~° 

(— 2 a — l)c 2 (a)ld(i) 

where B t is in G, a is the synchrotron spectral index 2 , K 0 is 
the ratio of the number densities of protons and electrons, I„ is 
the synchrotron intensity at frequency v, E p is the proton rest 
energy, and l is the path length through the radio relic. The 
constants ci, C 2 , and C 3 are defined as 

cr = 47 r ^ 3 e5 = 6.26428 x 10 18 erg -2 s _1 G _1 

C 2 = ^ 3 ( | 7 + + 7 y r[( 3 7 - 1 )/ 12 ] x r[(3 7 + 7)/12] 

/o 3 

c 3 = 77 - = 1-86 5 58 x 10 - 23 erg G _1 sr ” 1 

(47rm.eC 2 ) 

C 4 = [cos(i)]^ 7+1 ^ 2 , 



1/(3—a) 


,(4) 


where e is the elementary charge, m e is the electron mass, c is 
the speed of light, 7 is the spectral index of the electron energy 
spectrum which relates to the synchrotron spectral index a = 
— ( 7 — l)/2, and i is the inclination of the magnetic fields with 
respect to the sky plane (Beck & Krause 2005). 

We obtained averaged synchrotron intensity /„ of 6.27 x 
10 -19 erg s _1 cm -2 Hz -1 sr _1 at 2051 MHz in the radio relic 
where the fractional polarization was measured (the dashed 
frame region in Figure 1). On the other hand, we adopted the 
spectral index of a = —0.81 measured by van Weeren et al. 
( 2012 ) since we could not measure the accurate spectral index 
from our JVLA data due to the missing flux (see Section 4.2). 
We assumed the ratio of the proton-electron number densities 
of Ko = 100 , which is consistent with the acceleration process 
of cosmic-ray electrons for secondary particles and turbulence. 
Since the radio relic has a 25 kpc thickness at the minimum 
(Owen et al. 2014) and covers ~ 1125 kpc x 520 kpc (Clarke & 
Ensslin 2006), we assumed the path length l of 25 kpc and 1125 
kpc. For the inclination, we assumed a mid-value of i — 45°. 

We obtained the total equipartition magnetic field strengths 
of the radio relic of ~ 5.0 p,G with l = 25 kpc and ~ 1.8 pG with 
l = 1125 kpc. These values of micro-Gauss order are consistent 
with the values of 1.5lg ' 9 pG and 3 . 3 l 2 '2 pG with a = —1.25 
using the classical and the hadronic minimum energy condi¬ 
tions, estimated by Clarke & Ensslin (2006), respectively. 

To obtain the uniform and random magnetic field strengths, 
we can use the relationship between the observed fractional po¬ 
larization p and the degree of uniformity / of the magnetic fields 
(Segalovitz et al. 1976) 

3 7 + 3 [~ (l-/)7r 1/2 r[( 7 + 5)/4] r 1 

P 3 7 + 7 2/(sin6»)G+i)/2r[( 7 + 7)/4] ’ 

where T is the gamma function and 9 is the angle between the 


2 We use the definition of S oc v a . 
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line of sight and the uniform magnetic fields. The ratio between 
the strengths of the uniform magnetic fields B u and total mag¬ 
netic fields B t is given by B u /B t = f 2 ^ 1 + 1 ) (Beck 1982). The 
random magnetic field strengths B r is B r = (B 2 — B 2 ) 1 ' 2 . 

In order to avoid the effect of the depolarization, we used the 
fractional polarization of high frequency, which is p = 0.36 at 
3563 MHz. We also assumed the mid-value of 9 = 45°. 

We obtained the degree of uniformity of / ~ 0.56, which 
indicates that there are uniform magnetic fields in the radio 
relic with random magnetic fields. The uniform magnetic field 
strengths is ~ 3.7 pG with l = 25 kpc and ~ 1.3 pG with l = 
1125 kpc, and the random magnetic field strengths is ~ 3.4 pG 
with l = 25 kpc and ~ 1.2 pG with l = 1125 kpc. However, 
/ ~ 0.56 could be more large value because equation (5) does 
not take into account the depolarization. We can see the or¬ 
dered intrinsic magnetic fields over the beam size of ~ 52 kpc 
in Figure 6 against / ~ 0.56. This could indicate that there are 
random magnetic fields along the line of sight toward the radio 
relic and depolarization occurred. 


5.2 Magnetic Field Strengths in the Intracluster 
Space 


We estimated magnetic-field strengths in Abell 2256 using a 
traditional cell model (Lawler & Dennison 1982; Tribble 1991). 
In the model, we consider cells along the line of sight from 
the observer to the polarized source, and each cell consists of 
uniform size, uniform electron density, and uniform magnetic 
field strength with a single scale and random field orientation. 
In this case, distribution of RM becomes the Gaussian with zero 
mean, and the variance of RM is given by 

2 812 2 A b f, D ,2j, ,^ 

^- / {n e B ||) dl, (6) 


where Ab is the cell size in kpc. For the distribution of thermal 
electron density, we adopt the /3-model: 


n e = n 0 (l + -2 ) 


-3/3/2 


(7) 


where no is the central electron density in cm -3 , r is the dis¬ 
tance from the X-ray center in kpc, and r c is the core radius in 
kpc. Adopting equation (7) into equation (6), we obtain 


KBnorV 2 ^ 2 /r(3/3 — 0.5) 
°" RM_ (l +7 -2/ r 2)(6<3-:L)/ 4 V F(3/3) 


where B = y^By considering isotropic fields and V is the 
Gamma function. K is the constant which depends on the po¬ 
sition of a backside polarized source along the line of sight; 
K = 624 if the source is located behind the cluster and K = 441 
if the source is located at a halfway of the cluster (Feretti et al. 
1995; Felten 1996; Govoni et al. 2010). Therefore, with no, 
r, r c , (3 and Ab, equation (8) leads the magnetic field strength 
along the line of sight. 



Fig. 9. Histogram of the RM. We obtained 355 pixels in the RM map of 
Abell 2256. Dashed line represents the result of Gaussian fitting with 

/1 — —36 rad m — 2 and a = 7.4 rad m —2 . 


The adopted parameters and results are shown in Table 4. 
Here, we consider A s from 5 to 20 kpc according to a dynamo 
theory (Cho & Ryu 2009). We obtained that the field strength 
toward Source A is B = 1.26 pG with As = 5 kpc and B = 0.63 
pG with Ab = 20 kpc, and the field strength toward Source B is 
B = 0.21 pG with Ab = 5 kpc and B = 0.11 pG with As = 20 
kpc. 

5.3 Contribution of the Galactic Magnetic Fields 

Table 3 suggests a shift of the mean of the RM from 0 rad m -2 
to about —30 rad m -2 toward the Abell 2256 field. Figure 9 
shows the histogram of RMs in the Abell 2256 field. We ob¬ 
tained 355 pixels in Figure 6, and they actually indicates a his¬ 
togram centered around —36 rad m -2 . 

We consider that the shift is due to the Galactic contribu¬ 
tion. In order to estimate the Galactic contribution to the Abell 
2256 field, we examined the RM values of 28 polarized sources 
within 6° around Abell 2256 using the all-sky RM catalogue 
(Taylor et al. 2009), and calculated the average of their RMs. 
We found that the average of the RMs for 28 polarized sources 
is —30.0 rad m -2 with the standard deviation of 16.7 rad m -2 . 
The average is broadly consistent with the means of RM for the 
radio relic, Source A, and Source B. 

5.4 Step-like Variations of the Fractional Polarization 

The fractional polarization of the radio relic varies at ~ 3 GHz 
and around ~ 0.4-1.3 GHz, and gives step-like variations as 
shown in Figure 4. Such a decrease of the fractional polarization 
toward low frequencies implies that depolarization takes place. 
If this is the case, we could investigate turbulent magnetic fields 
along the line of sight, as introduced in Section 1. We can an¬ 
alytically model the fractional polarization in the cases of the 
EFD and IFD using the Burn’s law (Burn 1966). The fractional 
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Table 4. Parameters for magnetic field strengths in Abell 2256. 


Source 

K 

O’RM 

no* 

r t 

r c * 

P* 

Ab 

B 



[rad m -2 ] 

[10“ 3 cm- 3 ] 

[kpc] 

[kpc] 


[kpc] 

[fiG] 

Abell 2256 A 

441 

65.5 

2.6 

7.2 

587 

0.914 

20-5 

0.63-1.26 

Abell 2256 B 

441 

10.5 

2.6 

133.7 

587 

0.914 

20-5 

0.11-0.21 


* Reference: *Chen et at. (2007) ; ; Hireling et at. (1998). 


polarizations of EFD and IFD can be written as 

Pefd =Poe s , (9) 

and 

Pifd =Po — 77 —, (10) 

respectively, where po is the intrinsic fractional polarization, 
S = 2c>r M A 4 , (trm is the standard deviation of RM within the 
held of consideration, and A is the wavelength. The Burn’s law 
is, however, a function which does not produce a step-like vari¬ 
ation of the fractional polarization. Figure 4 shows an example 
of the EFD, and clearly indicates that a single depolarization 
component is not enough to explain the observed step-like vari¬ 
ations of the radio relic. 

In addition to the depolarization, we suspect that the missing 
flux also affects the fractional polarization above ~ 3 GHz. If 
there is large diffuse source larger than the LAS, we only detect 
the compact diffuse source smaller than the LAS. If the frac¬ 
tional polarization of the compact diffuse source is larger than 
that of the large diffuse source, then the fractional polarization 
could increase as the observing frequency increases. Therefore, 
the variation of the fractional polarization at 3-3.5 GHz may be 
partly due to the missing flux. 

Since the fractional polarizations of Sources A and B, which 
are located near the radio relic, does not show the variation at 
~ 3 GHz, its origin should be significant for the emission from 
the radio relic. Indeed, because Sources A and B are compact 
sources, the effect of the missing flux is expected to be insignif¬ 
icant for the emissions from Sources A and B. Thus, to clar¬ 
ify the effect of the missing flux on the fractional polarization, 
single dish observations at 3-3.5 GHz should be performed in 
future. 

On the other hand, the variation of the fractional polariza¬ 
tion expected at around ~ 1 GHz could not be explained by the 
missing flux, since the LAS at ~ 1 GHz or less is sufficiently 
larger than the the major axis of the radio relic. 

5.5 Depolarization toward the Radio Relic 

Hereafter, although the effect of the missing flux could be sig¬ 
nificant, we do not exclude the data above ~ 3 GHz in our anal¬ 
yses. This aims at studying the case that the variation of the 
fractional polarization at ~ 3 GHz is real. Again, we suggest to 
perform single dish observations at 3-3.5 GHz in the future, to 


clarify the effect of the missing flux on the fractional polariza¬ 
tion. 

Figure 9 implies that the histogram of RM follows the 
Gaussian distribution, which is not inconsistent to consider that 
the beam depolarization is induced by random (turbulent) mag¬ 
netic fields (Lawler & Dennison 1982). As already described in 
Section 5.4 and Figure 4, the Burn’s law with a single depolar¬ 
ization component is hard to reproduce the observed fractional 
polarization of the radio relic. Therefore, we consider models 
with multiple depolarization components along the line of sight 
toward the radio relic. A weakness of adopting the Burn’s law 
is that we cannot extract the information of magnetic fields. In 
order to understand the nature of depolarization as well as mag¬ 
netic fields, we hence carried out simulations of depolarization 
using simple random-field models. 

We calculated the polarized emission which pass the de¬ 
polarization components in the model. The components con¬ 
sist of grid, and each grid at the three dimensional coordinate 
(X,Y,Z) = (Nx Ny Nz) has single-scale random magnetic 
fields with uniform strength. The electron density is uniform 
in the components. The polarized intensity in each cylinder at 
(x,y) is given by 

P(x,y) = J poee^^’^dz (11) 

where po is the intrinsic fractional polarization and t is the syn¬ 
chrotron emissivity at a depth along the line of sight (Bum 
1966; Gardner & Whiteoak 1966; Sokoloff et al. 1998). The 
polarization angle rj>(x, y) is given by 

N Z 

(j>(x,y) = (j>o + 812^n e -B||A(A 2 (12) 

Z = 1 

where 4 >o is the intrinsic polarization angle in radians, n e is ther¬ 
mal electron density in cm -3 , By is magnetic field strengths 
parallel to the line of sight in pG, A l is the size of cells in kpc 
and A is observation wavelength in m. We obtain the polarized 
intensity P through the Nx x Ny cylinders as 

JV X Ny 

P = Y J Y j P{x,y). (13) 

x=\ y =1 

Our depolarization models include the following parameters: 
magnetic field strengths of a cell B, electron density n e , the size 
of the cells A l, the numbers of cells in the directions of the X, 
Y and Z axes (Nx, Ny, Nz), and the intensity of the polarized 
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Table 5. Parameters for the Depolarization Models. 


Model 

Component 

B 

[/+G] 

n e 

[1CT 3 cm -3 ] 

Al 

[kpc] 

Nx x Ny 
[kpc x kpc] 

N z 

[kpc] 

Intensity 

CRM 

[rad m —2 ] 

EFD+EFD 

foreside 

0.2 

1.0 

1* 

50 x 50 * 

600 

i 

3.2 


backside 

2.3 

3.0 



600 

3 

111 

IFD+IFD 

foreside 

0.2 

1.0 

5 * 

50 x 50 * 

500 

1 

4.1 


backside 

2.3 

3.0 



600 

3.5 

164 


* We assume 50 X 50 kpc since the beam size of 47 ;/ corresponds to ~52 kpc. 


(a) 


& 


EFD+EFD 



Fig. 10. (a) The depolarization model of EFD+EFD, which has non-emitting 
two depolarization components with two polarized sources, (b) The depo¬ 
larization model of IFD+IFD, which has emitting two depolarization compo¬ 
nents. In each model, the cells constituting the depolarization components 
have single-scale random magnetic fields with uniform strength and uniform 
electron densities. 



Frequency [GHz] 


source or emitting depolarization source. Adopting suitable pa¬ 
rameters, we can control the optimum frequency (Arshakian & 
Beck 2011) and the fractional polarization. 

We consider two components, EFD and IFD, and develop 
two two-component depolarization models, EFD+EFD and 
IFD+IFD (Figure 10). The order of the components along the 
line of sight from the observer is: 

• In the EFD+EFD model, we allocate the components in the 
order of a depolarization component, a polarized source, a 
depolarization component, and a polarized source. 

• In the IFD+IFD model, we allocate the components in the 
order of an emitting depolarization component, and another 
emitting depolarization component. 

Figure 11 shows the best-fits for the two models. Each pa¬ 
rameter is listed in Table 5. The EFD+EFD and IFD+IFD mod¬ 
els can nicely reproduce the fractional polarization of the radio 
relic. We confirmed that we need two depolarization compo¬ 
nents to produce the step-like variations of the fractional po¬ 
larization. In addition, we found that the <trm of the foreside 
depolarization component viewed from the observer has to be 
smaller than the that of the backside depolarization component. 

We could interpret the two depolarization components of 
the models as follows. The foreside depolarization component 
viewed from the observer could be the magneto-ionic plasma in 
the cluster or the Galaxy. The backside depolarization compo¬ 
nent should be magneto-ionic plasma inside the radio relic, ac¬ 
cording to the fact that the fractional polarization of Sources A 
and B does not show the same variation as the fractional polar¬ 
ization of the radio relic. Otherwise, if the cluster or the Galaxy 


Fig. 11. The fractional polarization spectra of the depolarization models. 
The solid and dotted lines show the EFD+EFD model and IFD+IFD model, 
respectively. Open squares represent the observed fractional polarization of 
the radio relic (the same as Figure 4). 

depolarizes relic’s polarization, the fractional polarizations of 
Sources A and B should also show the step-like variations. But 
as seen in Figure 5, the fractional polarizations of Sources A 
and B are nearly constant. The urm for the radio relic is larger 
than that for the cluster or the Galaxy. 

Although we can explain observed fractional polarization of 
the radio relic with our depolarization models as presented in 
this section, we should note that there should be more realistic 
model to match the observational data. For instance, Murgia et 
al. (2004) shows a realistic model obtained by numerical simu¬ 
lations, including the magnetic field strength, radial profile, and 
magnetic field power spectrum of cluster of galaxies, which is 
successfully applied to individual clusters (e.g. Govoni et al. 
2006; Guidetti et al. 2008; Bonafede et al. 2010; Vacca et al. 
2012 ). 

5.6 Faraday Tomography 

We also carried out Faraday tomography to make sure the exis¬ 
tence of multiple components toward the radio relic. We apply 
the so-called QU-fit in which a model is fitted with the data in 
Stokes Q and U spaces (see e.g., Ideguchi et al. 2014). We per¬ 
form the QU-fit using a Markov Chain Monte Carlo (MCMC) 
approach, so as to explore the best-set of model parameters. As 
for structures of Faraday components, we consider a delta func- 
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Fig. 12. Results of the QU-fit. Marks with error bars show the spatially- 
averaged fractional polarization data (Figures 4 and 11). Gray lines show 
the one-component fits (dashed: delta function, dotted: Gaussian). Black 
lines show the two-component fits (dashed: two delta functions, dotted: two 
Gaussians, solid: one delta function plus one Gaussian). 

tion or a Gaussian. Here, the delta function consists of the three 
parameters: the Faraday depth <j>, the amplitude, and the intrin¬ 
sic polarization angle xo, and the Gaussian function consists of 
the above three parameters plus the width of the Gaussian (the 
standard deviation of the normal distribution). We consider five 
models, one delta function, one Gaussian, two delta functions, 
two Gaussians, and one delta function plus one Gaussian, then 
find the best model according to the Bayesian information crite¬ 
rion (BIC). We also check the reduced chi-square (RCS) of the 
best fit for each model. 

The results are shown in Figure 12 and Table 6. We find that 
one-component models poorly reproduce the observed Q and 
U, and apparently the two-component models better fit with the 
data. Actually, two-component models dramatically improve 
BICs and RCSs (Table 6). We do not conclude the best-model, 
since the two-component models show similar BICs with each 
other. RCSs of the two-component models are a bit far from 
unity. To improve the fit, data below ~ 1500 MHz is essential. 

The fits with the two-component models commonly sug¬ 
gest that there are components at the Faraday depth <j> ~ 
—50 rad m -2 . This would be the radio relic, since the depth 
and the thickness are respectively close to the average and the 
standard deviation of RM for the radio relic (Table 3). 

6 Conclusions 

We reported new polarimetry results of Abell 2256 with JVLA 
at S-band (2051-3947 MHz) and X-band (8051-9947 MHz) in 
the C array configuration. We made images of the Stoke I, Q, 


and U, with 47" and 15. /r 1 resolutions. At S-band, we detected 
the significant polarized emission from the radio relic, Source 

A, and Source B. At X-band, we detected the significant polar¬ 
ized emission only from Source A. 

The total 2051 MHz flux of the radio relic is 286 mjy which 
includes the flux from the tail of Source C. The total flux is sub¬ 
stantially smaller than ~ 330 mjy, an expectation from previous 
L-band observations with an assumption of the spectral index 
a = —0.81. The estimated upper limit of the total 2051 MHz 
flux of the radio halo is ~ 49 mjy assuming that the radius of 
the radio halo emission is ~ 6'. 1 and the upper limit of the flux 
density is 606 pjy beam -1 . The estimated flux is also smaller 
than ~ 66 mjy, an expectation from previous L-band observa¬ 
tions with an assumption of the spectral index a = —1.1. 

We examined the missing flux caused by the largest angular 
scale (LAS) of our observations. Actually, the scale of the major 
axis of the radio relic and halo are 1014" and 732" (Clarke & 
Ensslin 2006) on the sky, respectively, while the LAS is 970" at 

1.5 GHz or 490" at 3.0 GHz in the JVLA C array configuration. 

We obtained RMs of the radio relic, Source A, and Source 

B. The mean and standard deviation of RM are (RM re ii c ) = 
—34.5 rad m -2 and <7RM reIic = 6.2 rad m -2 in the radio 
relic, (RMa) = —24.9 rad m -2 and ctrm a = 65.5 rad m -2 
in the source A, and (RMb) = —34.1 rad m -2 and <trm b = 

10.5 rad m -2 in the source B. 

We calculated the magnetic field strengths in the radio relic 
and the intracluster space in Abell 2256. For the radio relic, 
we calculated the magnetic field strengths using the revised 
equipartition formula, and the fractional polarization. The to¬ 
tal magnetic field strengths is ~ 5.0 pG with l = 25 kpc and ~ 
1.8 pG with l = 1125kpc, the uniform magnetic field strength is 
~ 3.7 pG with l = 25 kpc and ~ 1.3 pG with l = 1125 kpc, and 
the random magnetic field strength is ~ 3.4 pG with / = 25 kpc 
and ~ 1.2 pG with l = 1125 kpc. For the intracluster space, 
we calculated the magnetic field strengths using <trm- The 
magnetic-field strengths along the line of sight toward Source 
A is B ~ 1.26 pG with Ar = 5 kpc and B ~ 0.63 pG with 
Ar = 20 kpc. The magnetic field strengths along the line of 
sight toward Source B is B ~ 0.21 pG with Ab = 5 kpc and 
B ~ 0.11 pG with Ab = 20 kpc. 

We inferred that the shift of the mean of the RM from 
0 rad m -2 to about —30 rad m -2 toward the Abell 2256 field 
is due to the Galactic contribution since the averaged RM values 
of 28 polarized sources within 6° around Abell 2256 is broadly 
consistent with the means of RM for the radio relic. Source A, 
and Source B. 

We found that the fractional polarization of the radio relic 
remains about 20 % between 1.3-3 GHz and increases above 3 
GHz. The Burn’s law with a single depolarization component 
cannot reproduce the observed step-like fractional polarization 
spectrum of the radio relic. This may be due to the missing flux 
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Table 6. The reduced chi-square (RCS), the Bayesian information criterion (BIC), and best-fit values and 
1 -cr confidence regions for model parameters in the QU-fit. 


Model 

RCS 

BIC 


Amp. 

Xo 

Width 

Delta function 

21.2 

645.2 

—41.55^o 4 742 

0 49 OOO 8 

U.^tZ—0.008 

_n ££ 0.012 

U.OO_ 0.012 


Gaussian 

21.2 

648.6 

—41.53^ 3 7 6 75 

n 490008 

U.4Z_0.008 

n ££ 0 -oi 2 
— U.OD_o .011 

n 090 736 
U.UZ_o.o05 

two Deltas 

3.0 

110.2 

—40.87^ 2 o 6 16 

21 2 S —1,118 

-| 9 —0.003 
u ' iZ - 0.010 





-40.36^ 4 o 9 08 

21 17 -1,149 

Zl.l /_1 472 

— 1 45 -0 003 

i-^o_o.oio 


two Gaussian 

2.2 

93.3 

-49.32^i 9 0 1 03 

13.31^1 

0 49-0 013 
^^-0.051 

0 ofil-727 
O.OD_o.885 




-48.14^ 3 o 4 32 

13.4lJ;9Ji 

1 1 r— 0.012 

l-lO_ 0 .050 

0 7 rjl-913 
o-' U-0.733 

Delta + Gaussian 

2.6 

100.9 

-50.35°4 2 ? 6 2 

12 11 _1,270 
1 Z> 11 -1.555 

n ££0-069 
u - oo 0.045 





-49.38^f 2 0 ° 66 

1 9 9 Q — 1-244 
lZ.Zo _ 1 536 

1 noO-068 

— I.UZq.045 

1 £n 0 - 217 

1.0U_i.34i 


and/or depolarization. 

Our simulations of depolarization, which allow us to know 
three dimensional position information of the magneto-ionic 
plasma from the fractional polarization, indicated that two- 
component depolarization models can explain the step-like vari¬ 
ations of the fractional polarization. Furthermore, we found 
that the standard deviation of RM for the foreside component 
viewed from the observer should be smaller than that of the 
backside component. The existence of two components was 
also suggested from Faraday tomography. 

This work was supported by Japan Society for Promotion of 
Science KAKENHI Grant Numbers 26800104 (HN), 26400218 
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